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ABSTRACT 

The topology of the Las Campanas Redshift Survey is analyzed using Minkowski 
functional statistics, taking into account the selection effects of the survey. The results 
are compared with the predictions of some toy models of structure formation, including 
the Standard Cold Dark Matter Model and topological defect-based models. All of 
the toy models have a scale invariant primordial spectrum of perturbations, but quite 
different topologies. The statistics can discriminate between the predictions of the 
models with high significance. Amongst the four Minkowski functionals, the integrated 
mean curvature statistic appears to be the most powerful discriminant, followed by 
the genus statistic. None of the models considered gives an acceptable fit to the data. 
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At the present time the topology of the large-scale structure 
(LSS) of the Universe is not well known. Some surveys show 
pronounced coherent planar structures such as the "Great 
Wall" (see e.g. de Lapparent et al. 1986), others show evi- 
dence of filaments (see e.g. Haynes and Giovanelli 1986), and 
yet others appear dominated by superclusters and voids (see 
e.g. Kirshner et al. 1981). In many cases, however, these con- 
clusions are reached by visual inspection rather than by a 
rigorous statistical analysis. 

The topology of the LSS is very important for cosmol- 
ogy. Simple inflation-based models predict a Gaussian dis- 
tribution of the density perturbations (see e.g. Mukhanov et 
al. 1992, for a comprehensive review), whereas topological 
defect-based models predicted a non-Gaussian density field 
(see e.g. Brandenberger 1994, for a review). These models to 
a first approximation all predict a scale-invariant spectrum 
of primordial perturbations; and therefore the predictions of 
the models cannot be differentiated by means of power spec- 
tra and two-point functions. However, defect models can be 
differentiated amongst one-another as well as from inflation- 
based models in terms of the topology of LSS which they 
yield. 

The currently favored simple inflation-based theories 
predict a Gaussian distribution of the density fluctuations. 
Nonlinear evolution and biasing are expected to produce 
secondary non-Gaussianities, but of a type very different 



from the primary non-Gaussian signals predicted in defect 
models. For example, the genus curve (one of the Minkowski 
functionals used in our work) of a simulated SLOAN redshift 
survey is indistinguishable from that of a Gaussian distribu- 
tion (CoUey et al. 1999). 

The nontrivial topological signatures of most defect 
models are impossible to detect in the current CMB exper- 
iments (and hard to detect even in future experiments such 
as MAP) since due to the central limit theorem the non- 
Gaussian effects disappear at the large angular scales ac- 
cessible (see e.g. Moessner et al. 1994, for a specific study). 
Conversely, the topology of the LSS may be less sensitive 
to complicated biasing effects than the topology on smaller 
scales, thus rendering the study of the topology of the LSS 
an interesting endeavor. 

Recently, the use of the Minkowski functional statistics 
has been suggested as a useful way to study the topology 
of structure in the Universe (Mecke et al. 1994, Schmalz- 
ing and Buchert 1997). It has been shown (Mitsouras 1998, 
Mitsouras et al. 1999) that Minkowski functionals can dif- 
ferentiate between models with identical power spectra such 
as defect toy models and infiation-based models. Amongst 
the four Minkowski functionals MO — M3, the M2 statis- 
tic (integrated mean curvature) was shown to be a better 
discriminator than the frequently used genus statistic M3. 
Minkowski functionals have been applied to Abell galaxy 
clusters surveys (Kerscher et al. 1997) and to the IRAS cat- 
alog (Kerscher et al. 1998), among others, and to numerical 
simulations (Schmalzing et al. 1999). 

Another set of useful statistics to explore the mor- 
phology of the large-scale structure are the structure func- 
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tions (Babul and Starkman 1992, Sathyaprakash et al. 1996) 
which have e.g. been applied to the 1.2 Jy IRAS catalog 
(Sathyaprakash et al. 1998). These statistics are based on 
ratios of eigenvalues of the inertia tensor of the mass dis- 
tribution and carry information different from that of the 
Minkowski functionals. On the other hand, the Shapefinder 
statistics (Sahni et al. 1998) are ratios of Minkowski func- 
tionals. They have very recently (Bharadwaj et al. 2000) 
been used to analyze the LCRS data set (without taking into 
account the important selection effects), finding evidence for 
filamentarity. 

The deepest wide field redshift survey to date is the 
Las Campanas Redshift Survey (LCRS) (Shectman et al 
1996). An initial analysis (Doroshkevich et al 1996) indi- 
cated that the LCRS field was dominated by planar struc- 
tures, with higher density filaments at their intersections. 
In contrast, the analysis of Mitsouras (Mitsouras 1998) us- 
ing the Minkowski functional statistics detected no devia- 
tions from a Gaussian distribution. Neither analysis, how- 
ever, took into account the complicated selection effects in- 
trinsic to the LCRS survey. 

In this paper, we calculate the Minkowski functionals 
for the LCRS, taking careful account of the nontrivial se- 
lection effects of that survey (except for the limited surface 
brightness selection effect). We compare the resulting statis- 
tics with the predictions of some toy models of structure for- 
mation: the Standard Cold Dark Matter (SCDM) model as 
a representative of the class of inflation-based models, and 
three toy defect models (based on string wakes, string fila- 
ments, and textures, respectively). Our aim in this paper is 
to focus on the potential of the LCRS survey to differentiate 
between different models of structure formation with iden- 
tical power spectra rather than on studying specific COBE 
normalized models. We find that no model comes close to 
the data. In work in progress we are studying a class of 
COBE-normalized inflation-based models with the goal of 
discovering whether some members of this class come closer 
to explaining the observed topology than the SCDM model. 



2 DATA ANALYSIS 

The Las Campanas Redshift Survey (LCRS) is a fiber-optic 
survey consisting of 26418 redshifts of galaxies. 1.5° x 1.5° 
spectroscopic fields were observed to pave 700 deg^ of sky 
in six strips, each roughly 1.5° x 80°. Three strips centered 
on 5 = —3°, —6°, and —12° are in the North galactic cap 
(NGC) and three centered on 5 ^ -39°, -42°, and -45° 
are in the South galactic cap (SGC). The LCRS probed to 
redshifts z ~ 0.2 with mean redshift z ~ 0.1 or equivalently 
to distances of 600 /i~^Mpc. Approximately 70% of galax- 
ies having magnitudes in the range 15.0 ^ rrir < 17.7 were 
measured spectroscopically. For our analysis, we only con- 
sider galaxies out to a redshift z < 0.15 or distances of 450 
h~^Mpc or less. 

The selection criteria for the LCRS are more complex 
than a simple magnitude limit and have their origin in the 
technique of the fiber-optic measurements. Galaxies were 
selected from a CCD-based catalog based on four selection 
criteria. For the first 20% of survey galaxies, redshifts were 
measured using a 50-fiber spectrograph system and the fi- 
nal 80% using an improved 112-fiber system. Galaxies were 



first selected by applying both faint and bright isophotal 
magnitude limits, with nominal values 

mi ^ m < 7712 : 

r mi = 16.0, m2 = 17.3, 50 - fiber data (1) 
I mi = 15.0, m2 = 17.7, 112 - fiber data ' 

The second criterion selects against low surface brightness 
(LSB) galaxies by imposing the nominal central surface 
brightness cutoff 

rric < rricen — 0.5(m2 — m) : 

J nicen = 18.15, 50 — fiber data (2) 
1 rricen = 18.85, 112 — fiber data ' 

where the central magnitude mc corresponds to the flux 
within a 2 pixel radius of the centroid of the galaxy. There 
are small field-to-field variations from the nominal photo- 
metric limits which were made in response to variations in 
sky brightness and in seeing. The third criterion restricts the 
maximum number of galaxies in a given field to the num- 
ber of fibers available during observation of that particular 
field. Lastly, the fourth criterion imposes a minimum separa- 
tion of 55" between galaxies because mechanical constraints 
prevent the fibers from being placed any closer together. 



3 MODELS 

There are two main classes of theories for structure forma- 
tion which may possibly be able to account for the origin of 
primordial perturbations and observed large-scale structure 
in the universe. The first are infiation-based models in which 
quantum fluctuations are generated at an early stage of in- 
flation on microphysical scales, and the wavelength of these 
perturbations subsequently grows exponentially and thus 
gives rise to density perturbations on cosmological scales 
(Mukhanov and Chibisov 1981, see also Press 1980, Lukash 
1980, Chibisov & Mukhanov 1980, Sato 1981, Chibisov & 
Mukhanov 1982 for related work). Simple inflation-based 
models predict a scale-invariant spectrum of adiabatic per- 
turbations with random phases, each having a Gaussian dis- 
tribution. The second class is based on topological defects, 
formed during a symmetry breaking phase transition in the 
early universe, giving rise to seeds for structure formation 
(Kibble 1976, Zel'dovich 1980, Vilenkin 1981). Topological 
defects represent regions in space with trapped energy den- 
sity and they result in a roughly scale-invariant spectrum of 
non-adiabatic perturbations with non-Gaussian phases. In 
this paper, we compare the inflation-based "standard" cold 
dark matter (SCDM) model with topological defect models 
like cosmic strings and textures. ^ 

* After this work was completed, the results of the 
BOOMERANG (de Bernardis et al 2000) and MAXIMA-I 
(Hanany et al. 2000) CMB experiments were announced, which 
with high statistical significance demonstrate the existence of a 
well-defined and narrow first acoustic peak in the angular power 
spectrum, thus ruling out topological defects as the main source 
of density fiuctuations (see e.g. Albrecht et al 1999 and references 
therein). It is, however, important to keep in mind that defects 
may still contribute a significant amount to the primordial power 
spectrum (see e.g. Bouchet et al. 2000), and that the resulting 
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3.1 Standard Cold Dark Matter Model (SCDM) 



v^{ti)z{ti) X ^{ti)z{ti) X h{ti) 



(4) 



In the SCDM model, the dominant constituent of the uni- 
verse is cold dark matter and the standard model is do- 
scribed by a density parameter ft = 1, Hubble constant 
Ho = lOO/i km/s/Mpc with h = 0.5, and cosmological 
constant A = 0. This specific model confiicts with obser- 
vations if we compare the normalization of the power spec- 
trum of density perturbations with that of the CMB fiuctu- 
ations. Other inflationary models like tilted CDM, ACDM, 
and mixed dark matter models have the potential to ac- 
count for observations. Since we arc primarily concerned 
with the power of the Minkowski functionals to differenti- 
ate between the signatures of primordial Gaussian and non- 
Gaussian phases, we have chosen to work with the standard 
model. We have also modified the standard value of the Hub- 
ble constant from h = 0.5 to h = 0.7 to reflect some of the 
more recent measurements. 

In the SCDM simulation, 192^ particles are evolved 
in a cube of comoving length of 630 h~^Mpc using the 
Bertschinger & Gelb (1991) P'''M N-body code. The normal- 
ization of the primordial power spectrum was chosen such 
that CT8 = 0.95 to agree with the observed clustering (Pen 
1998). A COBE normalized model would require era = 1.22. 



3.2 Cosmic Strings 

Cosmic strings are linear topological defects which arise if 
the vacuum manifold Ai is not simply connected. At the 
time of the phase transition, a network of strings is formed. 

The phase transition takes place on a short time scale At < 
, and will lead to correlation regions of radius $, < t. The 
strings are either infinite in length or closed loops. Infinite 
strings give rise to filaments or wakes while loops form seeds 
for spherical accretion of matter and radiation. In our study, 
we are concerned only with the infinite strings, but realistic 
cosmic string models should contain a mixture of both. 

In the string wake model, strings are straight on the 
scale of the Hubble radius and the string tension |p| equals 
the mass per unit length /i. As a result, the strings will 
move relativistically and exert no local gravitational force. 
The strings will have a large transverse velocity v and give 
rise to a planar velocity perturbation. A string at time U 
produces a wake of comoving length v(,{ti)z{ti) and width 
^{ti)z{ti). The wake thickness can be calculated using the 
Zel'dovich approximation (Zel'dovich 1970) and will depend 
on the nature of the dark matter since free streaming has an 
effect. For hot dark matter, the thickness h{ti) at the present 
time to of a wake created at time ti > teq is (Perivolaropou- 
los et al. 1990, Brandenberger 1991) 



24n 



h{ti) = ^^-^Gtivy{v)z(tiY^^to 



(3) 



where 'y{v) is the relativistic factor associated with v. The 
cosmic string model of structure formation is viable even 
if dark matter is hot since free streaming will not erase the 
non-adiabatic primordial perturbations on a small scale. The 
wakes will have dimensions 



non-Gaussian signals in the large-scale structure of the Universe 
may well be important. 



1/3,8/3 



(5) 



with mass 

m{ti) = '^^Gfiv^^{v)u~^potl 

The string correlation length is taken to be ^{t) = t/v, 
where f could in principle be determined from numerical 
simulations. Hot dark matter wakes created before teq arc 
diluted by the radiation pressure and, therefore, do not con- 
tribute any significant amount to structure formation. The 
string waJie model predicts that planar wakes are the dom- 
inant structures in the large-scale galaxy distribution with 
the most rmmcrous and thickest wakes created at teq. 

In the string filament model, strings have a lot of small- 
scale structure and the string tension \p\ will be smaller than 
/i. The strings will have a small transverse velocity and will 
exert local gravitational attraction. This induces a filamen- 
tary pattern of accretion. Zanchin et al. (1996) have ana^ 
lyzed the string filament model for moderate transverse ve- 
locities and we will here present an extreme case in which the 
transverse velocity is negligible (Mitsouras et al. 1999). Us- 
ing the Zel'dovich approximation, the comoving radius q{ti) 
of a filament created at time ti is (Aguirre & Brandenberger 
1995) 



q{ti) = 



12 



GAln 



1/2 



to , 



(6) 



where A = // — The mass in an individual filament is then 

(7) 



U\ l-^TT _i ,1/3,8/3, /to 

m{ti) = ——GXp Pot/ V In — 

where v is defined £is previously. 

In the toy topological defect simulations, galaxies are 
just laid down and no actual evolution takes place. The pro- 
cedure for setting up a string filament or string wake model 
first required dividing the time interval between teq and to 
into Hubble times. At each time step ti, a network of strings 
with correlation length ^{ti) is laid down to simulate the evo- 
lution of the strings. For each time step ti, beginning with 
teq, we divided the simulation cube into N{ti) Hubble vol- 
umes ti - A fixed number cj^ of strings per Hubble volume is 
laid down, choosing string centers and directions at random 
in the cube. Thus at each time step, cj^ x N{ti) is the num- 
ber of strings being laid down in the cube. For given values 
of GA and i' in the case of the string filament simulation and 
Gfi, iJ and V in the case of the string wake simulation, the 
total number Cs of strings is determined by demanding that 
the total mass satisfies critical density. Galaxies are then laid 
down at random in the volume of the strings, with a uniform 
density amongst the strings. For our string filament simu- 
lation, we have taken GA = 10~® to agree with the COBE 
normalization and f = 1, such that the length of the string 
is equal to the Hubble radius at its time of creation. The 
string wake simulation has values G/tt = 10~^ and = 1. 



3.3 Textures 

Textures are topological defects whoso vacuum manifold M 
is a sphere . In 4-dimensional space-time, textures are 
unstable (Turok 1989) and can collapse to form seeds for 
structure formation. While the energy in cosmic strings is 
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trapped in potential energy, in textures, the energy is in 
spatial gradients with possibly some in kinetic energy. 

The simplest texture configurations are spherically sym- 
metric (Turok 1989). At each time ti, there is a finite proba- 
bility p that in a volume tf, there will be a texture configura- 
tion entering the horizon (Prokopec 1991, Leese & Prokopec 
1991). The initial texture configuration cam be realized by 
a spherically symmetric velocity perturbation which results 
in a spherically symmetric density perturbation (Gooding 
et al. 1991). The comovirig radius q{ti) of this perturbation, 
calculated (Aguirre 1995) using the Zel'dovich approxima- 
tion, is 



1/3 4/3 



(8) 



where e = Gr;^ and r] is the scale of symmetry breaking. 

In the toy texture model, p spherical texture configura- 
tions are laid down per Hubble volume tf with centers chosen 
at random in the box. Galaxies are laid down within each 
texture with a Gaussian radial density profile and random 
angular distribution. The number of galaxies per texture 
was chosen to satisfy critical density for the entire volume. 
A value of e = 3 x 10® is an appropriate value for a COBEl- 
normalized texture model (Pen et al. 1997). 



3.4 Construction of mock LCRS catalogs 

We construct mock catalogs by applying the LCRS selection 
criteria to the simulations. The first step involved randomly 
assigning luminosities and corresponding absolute magni- 
tudes to all the galaxies in each cube, according to the 
Schechter luminosity function 



-L/L' 



(9) 



where d>* = 0.019 h' Mpc ' , a = —0.70, and L* corresponds 
to the absolute magnitude M* = —20.29 -|- 5 log h for the 
LCRS (Lin et al. 1996). The faint end slope a = -0.70 
results from fitting the LCRS luminosity function and not 
from the absence of intrinsically faint galaxies in the sur- 
vey. The next step was to assign a surface brightness to the 
galaxies, but this presented a significant problem since the 
distribution of surface brightness is nontrivial. While LSB 
galaxies are found in the same large scale structure as other 
galaxies, they tend to be less strongly clustered (Mo et al. 
1994). Since biasing may occur in one or more of the models, 
the LSB selection criterion has been omitted. 

To apply the selection criteria, slices of appropriate ge- 
ometry are first cut from the simulation cubes. We have 
chosen to consider the three strips in the North galactic 
cap as a composite slice and the three strips in the South 
galactic cap as another composite slice. The 3-dimensional 
geometry of the NGC and SGC slices should reveal more 
information about the topology of the galaxy distribution. 
To avoid boundary problems in the 630 /i^^Mpc simulation 
cube, slices are cut from an interior cube with side-lengths of 
490 ft~^Mpc. After converting the absolute magnitudes into 
apparent magnitudes, the exact faint and bright apparent 
magnitude limits in each spectroscopic field (Sliectman ct 
al. 1996) are applied. Next, the restriction on the maximum 
number of galaxies in each spectroscopic field is imposed, 
depending on the fiber system used for that particular field 



(random elimination of galaxies if the number of galaxies in 
the simulation field exceeds the maximal number of fibers 
used for that field). Finally, to apply the minimum sepa- 
ration criterion without introducing bias, the galaxies are 
ordered in RA to form the set {gi\gi, giv}. Starting with 
galaxy gi and choosing to keep it, any galgixy within 55" is 
immediately removed from the set. This step is iteratively 
carried out until the end of the set is reached. 



4 MINKOWSKI FUNCTIONALS 

Hadwiger's theorem states that in d spatial dimensions, the 
global morphological description of any spatial distribution 
can be completely characterized by d+1 descriptors com- 
prising the set of Minkowski functionals (Hadwiger 1957). 
Given an isodensity surface dK in 3d, the morphology is 
completely determined by the following four statistics. The 
first Minkowski functional 

Mo = y , (10) 

is simply the volume enclosed by the surface dK and the 
second functional 

Ml = I , (11) 

is related to the surface area A. The third functional 

M. = A , (12) 

provides information about shape and is proportional to the 
integrated mean curvature 



(13) 



where Ri and R2 are the principal radii of the surface dK 
The fourth Minkowski functional 

3X 



M3 



47r ' 



which is related to the Euler characteristic 
1 



dK 



RiR'. 



-dA . 



(14) 



(15) 



is purely a topological quantity and can be geometrically 
defined as 

X oc (# of clusters) — (# of tunnels) -|- (# of cavities) .(16) 

The Minkowski functionals are realized for a galaxy dis- 
tribution by placing spheres of radius r centered on the 
galaxies. The union of these spheres specifics an ensemble 
of clusters, where two spheres are said to belong to the 
same cluster if they are connected by a chain of intersecting 
spheres. The Minkowski functionals will be used to measure 
the topology of the coverage. We use the method of Mecke 
ct al. (1994) and Kerscher et al. (1997) and display the func- 
tionals as a function of the length scale r, which acts as a 
smoothing length of the density field. 



5 RESULTS 

For each model studied, the simulations were performed in 
the following way. First, a cubical box of comoving length 
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630 /i-^Mpc was populated with galaxies according to the 
prescriptions discussed in Section 3. Next, slices with the 
appropriate geometry of the LCRS survey were cut out from 
the simulation cube. To avoid boundary effects, the slices 
were extracted from an interior cube with side length of 490 
/i~^Mpc. For comparing with the individual LCRS slices, 
four galaxies in the interior box were chosen as centers of 
slightly oversized slices of 90° in length and 1.6° in width. 
For comparing with the composite LCRS slices, four larger 
slices of 90° in length and 10.6° in width were cut out. 

Luminosities were assigned randomly in the initial sim- 
ulation box according to the Schechter luminosity function. 
For each slice, this then allows us to determine the apparent 
magnitude of all galaxies contained in the slice. Next, the 
LCRS selection criteria were applied as detailed in Section 
3.4. 

For each slice, the Minkowski functionals were calcu- 
lated using a slightly modified version (Mitsouras 1998) of 
the program package publically available (Kerscher et al. 
1998). Averages over the four centers were taken. Under the 
ergodic hypothesis, the volume averages (averages over the 
slices) will be equal to the ensemble averages over different 
realizations of the theories. The error bars indicated in our 
plots correspond to these statistical errors. 

Figures 1-4 show the resulting volume- averaged 
Minkowski functionals for the NGC slice. The results of our 
simulations are compared with the LCRS data. As men- 
tioned in the Introduction, none of the models considered 
in this paper is able to explain the data. The SCDM model 
gives the least bad fit, but even it is ruled out with high 
statistical significance. |^ 

As is apparent from Figures 1-4, the Minkowski func- 
tionals provide a very powerful tool to distinguish between 
models with identical power spectra but different topology. 
Amongst the Minkowski functionals, the integrated mean 
curvature statistic M2 appears to be the most powerful dis- 
criminant, followed by the genus statistic M3. This confirms 
the results seen in previous studies (Mitsouras 1998, Mit- 
souras et al. 1999, Kerscher et al. 1997). 

The predictions of the texture toy model differ dra- 
matically from those of the other models studied, and also 
from the data. The reason is that in the texture toy model, 
structure is dominated by well separated spherical clumps, 
whereas in the other models there is a higher degree of con- 
nectedness of the individual structures. In addition, the cen- 
tral density in the texture toy model is larger than in the 
other theories studied. Therefore, MO, the volume enclosed 
by the "isodensity" surface, is significantly smaller for any 
value of r than in the other models. The surface area (Ml) 
itself continues to increase with r since to a first approxi- 
mation the number of disconnected components remain the 
same while the size of each object increases. In the other 
models studied, the value of Ml for small r is higher since 
the shape of the objects does not minimize the surface area 
for fixed enclosed volume. However, as r increases, the grow- 
ing interconnectedness of the components leads, for large r, 
to a gradual decrease in the value of Ml. The radius corre- 



t The oscillations seen at large values of r are artifacts of the 
boundary conditions. 



sponding to the maximum of Ml corresponds to the mean 
separation between clusters. 

For a single component surface without holes and with 
fixed volume, the sphere minimizes the integrated mean cur- 
vature. Hence, for small values of r, the value of M2 is lower 
in the texture model than in the other models. In particular, 
in the string filament and wake toy models the individual 
filaments and wakes have a small curvature radius in one 
direction, and hence a large mean curvature. Once the indi- 
vidual filaments and wakes merge by overlapping, the overall 
shape of the surface becomes concave and hence M2 turns 
negative. This feature is shared by the SCDM model, and by 
the data. Similarly, the overlapping of filaments and wakes 
leads to a structure with many holes for larger values of r. 
This leads to negative genus for these scales, a feature which 
is again shared by the SCDM model and by the data, but is 
absent in the texture toy model. The radius corresponding 
to the minimum of M3 gives another good indicator of the 
mean separation of clusters. 

By comparing the LCRS data with the simulations, it 
appears that the data has less connectivity than what would 
be seen in the SCDM model, and significantly less than what 
would be predicted by the toy string theories. These results 
are confirmed by the analysis of the SGC slice (Figures 5 
- 8). Note that the differences between the NGC and SGC 
slices are due to the different selection effects. 

The number of galaxies in the model slices is not the 
same as in the LCRS slices. In order to check that the dis- 
agreement between the theories and the data is not due to 
this difference, we have repeated the analysis with "ran- 
domly sampled" model slices, model slices from which in 
a random way galaxies were removed in order to obtain the 
same number of galaxies as in the actual LCRS slices. Fig- 
ures 9-12 give the results for the "randomly sampled" NGC 
slice. Figures 13 - 16 the corresponding results for the SGC 
slice. As is apparent, the conclusions are unchanged. 

As mentioned in Section 3.4, we did not apply the LSB 
selection criterion of the observational survey. To investigate 
whether this omission could account for the disagreement 
between the data and the SCDM model, we repeated the 
calculation of the Minkowski functionals of the observational 
data with the lowest 20% of the galaxies (in terms of surface 
brightness) removed. There were no appreciable changes in 
the results. 



6 DISCUSSION 

The main conclusion of our study is that Minkowski func- 
tionals provide a very powerful statistical tool to differen- 
tiate between theories with identical power spectra but dif- 
ferent topologies. We have compared the predictions of the 
SCDM model and of three topological defect toy models 
with the LCRS data, taking careful account of most of the 
selection criteria of the LCRS survey. 

None of the models studied comes close to giving a sat- 
isfactory match to the data. This is not surprising since the 
SCDM model is known not to correctly describe the ob- 
served clustering of matter in the Universe, and since the 
defect models considered here are all unrealistic toy mod- 
els. The texture model produces isolated and well-separated 
dense clusters and gives a picture of the distribution of 
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galaxies which already at first glance is inconsistent with 
any observational data set. The string models are based on 
ideal straight filaments or wakes, and the filaments and walls 
seen in observations are not straight. 

It would be of great interest to explore whether more 
realistic models can provide a better fit to the data. In 
the class of infiation-based models yielding a spectrum of 
roughly scale-invariant primordial perturbations it would be 
of interest to explore the effects of reducing Q.m and intro- 
ducing a cosmological constant A. In particular, it would be 
of great interest to see if the Minkowski functionals as ap- 
plied to the LCRS geometry can differentiate between the 
set of infiation-based models which are compatible with the 
present CMB data and with the amplitude of the matter 
power spectrum on large scales. Work on this issue is in 
progress (Merali et al. 2000). 

From the point of view of particle cosmology, the defect 
models studied in this paper are too idealized to be realistic. 
In texture models the resulting structures are not all spher- 
ically symmetric, and there are important non-topological 
fiuctuations which also lead to density inhomogeneities (Bor- 
rill 1994). Similarly, the structure seeded in a cosmic string 
model will also include effects of string loops (see e.g. Bran- 
denberger 1994 for a review), and the structures seeded by 
the long string network will be given by a combination of 
the filament and wake effects. In addition, the structures 
will not be straight. It would also be of interest to use the 
Minkowski functionals to study the predictions of numeri- 
cally simulated more realistic defect models, and to see if 
there are parameters for which good agreement with the 
data is obtained. 

It is also possible that the differences between the theo- 
retical predictions and the LCRS data are due to biasing ef- 
fects. This issue could be investigated by using hydrodynam- 
ical simulations of structure formation. Redshift distortions 
also complicate the interpretation of clustering in redshift 
surveys. The peculiar velocities of galaxies cause them to ap- 
pear displaced along the line of sight in redshift space, thus 
causing the observed galaxy distribution to be anisotropic. It 
would be interesting to compare the real and redshift space 
results from N-body simulations using Minkowski function- 
als. 

In conclusion, we hope to have convinced the reader 
that Minkowski functionals provide a very powerful tool for 
quantitative cosmology. 
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Figure 1. MO NGC slice without random sampling. 




10 

r [h"^ Mpc] 



Figure 3. 



0.009 



0.008 



0.007 



0.006 



0.005 



0.004 



0.003 



0.002 



0.001 - 



10 

r [h"'' Mpc] 



LCRS 

CDM ^B- 
Filament — ^ 
Wake — ©- 
Texture — e- 




3e-05 



2e-05 



1e-05 



- 



-1e-05 



-2e-05 



-3e-05 




Figure 2. 



r[h"' Mpc] 
Figure 4. 



© 0000 RAS, lUNRAS 000, 000-000 



8 H. Trac et al. 




10 

r [h"^ Mpo] 

Figure 5. SGC slice without random sampling. 



0.0008 



0.0007 - 



0.0006 - 



0.0005 



0.0004 



0.0003 - 



0.0002 



0.0001 



1 

LCRS 
CDM 
Filament 
Wake 
Texture 



-B- 
-0- 

-e- 




Figure 7. 



0.009 



0.008 



0.007 



0.006 



0.005 



0.004 



0.003 



0.002 



0.001 - 



10 

r [h"'' Mpc] 



LCRS 

CDM ^B- 
Filament — ^ 
Wake — ©- 
Texture — e- 




4e-05 



3e-05 



2e-05 



1e-05 



-1e-05 



-2e-05 



-3e-05 



-4e-05 




Figure 6. 



r[h"' Mpc] 
Figure 8. 



© 0000 RAS, MNRAS 000, 000-000 



Minkowski Functionals 9 




r [h"^ Mpo] r [h"^ Mpc] 



Figure 9. NGC slice normalized with respect to the average num- Figure 11. 

ber of galaxies in the LCRS NGC slice. 




r [h"^ Mpc] 

Figure 12. 

Figure 10. 



© 0000 RAS, MNRAS 000, 000-000 



10 H. Trac et al. 




10 

r [h"^ Mpo] 

Figure 13. SGC slice normalized with respect to the average 
number of galaxies in the LCRS SGC slice. 
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